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Abstract 



The main purpose of this work is the investigation of the influence of clcctron- 
positron-ion on the electrostatic wakefields that are deriven by intense, short 
electromagnetic wave packets in a three component unmagnetized plasma. The 
equations that describes the deriven electrostatic wakefield are given in the form 
of system of nonlinear differential equations. The solution of that system was 
carried out numerically and the results are displayed graphically. 



Introduction 



Recently, the nonlinear propagation of electromagnetic waves in electron-positron 
plasma has attracted the interest of researchers Q [|j and|3|]. Due to the fact 
that the electron-positron plasma are found in Van Allen Belts, near polar cap 
of pulsar, in the active galactic nuclei, as well as in the early universe. 

When plasma becomes so hot that it becomes relativistic the temperature 
T of the plasma exceeds the rest mass energy of electrons mc 2 = 0.5 MeV. 
In this relativistic regime the processes of electron-positron pair creation and 
annihilation become important; 27 ^ e + + e~ . In relativistic temperature the 
electrons (and positrons) energy E e far exceeds the rest mass energy so that 
electrons and positrons behave kinematically similar to photons and come into 
equilibrium with nearly equal population. 

In 1990 T. Tajima and T. Taniuti pj are interested in the process of electron- 
positron plasma creation and annihilation occurs in relativistic; unmagnetizcd 
plasma at high temperature of the plasma exceeds the rest mass energy of 
the electrons. They are emphasize the dominant population of electrons and 
positrons and neglect ions effect. 

In || , it was shown that , the positrons can be used to probe the parti- 
cle transport in Tokamak plasma. This new diagnostic technique could also 
be useful in studying transport in other magnetic confinement devices such 
as reversed-field pinches and magnetic mirrors. Electron-positron pair pro- 
duction can also be possible during intense short laser pulse propagation in 
plasma. When picosecond electromagnetic pulse having intensity beyond 10 21 
W/cm 2 interact with plasma, the plasma electrons can acquire oscillatory ve- 
locity (V os = eE /mc) higher than the speed of light c. Here E (W ) is the 
amplitude (frequency) of the EM radiation, c & m are respectively the magni- 
tude of the electron charge and the rest mass of the electron. 

Several authors have suggested that when the oscillatory electron energy 
v 2 

E os = mc 2 (\ + -^f) 1 ' 2 exceeds 3mc 2 then these electrons can create an electron- 
positron pairs in the presence of background positive ions. And also it has been 
shown that powerful short laser pulsar can generate a large amplitude electro- 
static wave (Wakefield) in a plasma in which the electron plasma frequency u> pe is 
much smaller than io Q . Since the amplitude of the excited wakefield can be very 
high ( e4>/m c 2 = $ ^> 1, where <j) is the electrostatic potential ), the plasma 
electrons also acquire the longitudinal energy El = m c 2 (l + / (m c 2 )) ^> 
3m c 2 , where Pm is the relativistic momentum in the wakefield. Thus, there 
is the possibility of creating electron-positron pairs in those regions in which 
wakcfields and short EM pulses are localized. Since the lifetime r p > lj^ 1 of 
the positron in such a plasma is sufficiently long. 

We shall have a plasma that is an admixture of electrons, positrons and 
positive ions. A three-component electrons-positron-ion plasma can indeed be 
created in laboratory plasma. V. I. Berezhiani, and et. al. || have investigated 
the influence of electron-positron pairs on the electrostatic wakefield that are 
driven by intense, short EM wave packets in a three component unmagnetizcd 
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plasma. They found that in contrast to pure electron-ion plasma, the pres- 
ence of positrons significantly reduce both the amplitude and the wavelength of 
wakefields. In view of their investigation, it is suggestive that the production 
of electron-positron pairs by wakefields in Laboratory and cosmic plasma ought 
to be reconsidered because the presence of pairs may affect the wakefield that 
is responsible for accelerating electrons to very high energy. 

The important of three component admixture plasma had led to much re- 
lated theoretical investigations. Of these, Rizzato @, investigated the possibility 
of high-frequency EM localization in cold unmagnctized plasmas made up of 
electrons, positrons and ions. It was shown that such a possibility depends on 
the concentration of particles, on the velocity of the solitons and on the angle 
between the direction of modulation and the direction of the fast spatial de- 
pendence. In particular in the case of pure electron-positron plasmas Rizzato 
was shown that no localized solution is possible which is in contrast with results 
derived by Mofiz et. al. (1985) from a similar model. Also it was shown that 
low frequency magnetic fields may appear in the case of oblique modulation 
and that both the amplitude of the magnetic field and the amplitude of the 
solitary wave are very sensitive functions of the angle between the direction of 
modulation and fast spatial dependence. 

Berezhiani, et. al. Q analyzed the nonlinear interaction of an arbitrarily 
large amplitude circularly polarized EM wave with an unmagnctized electron- 
positron plasma, taking into account relativistic particle-mass variation as well 
as large-scale density perturbation created by radiation pressure. It was found 
that the interaction is governed by an equation for the electromagnetic wave 
envelope, which is coupled with a pairs of equations describing fully nonlinear 
longitudinal plasma motion. This results should be useful for the understanding 
of nonlinear photon motion in cosmic plasmas, such as those found in the early 
universe and active galactic nuclei. 

Mofiz 0] studied the nonlinear propagation of Langmiur wave in a hot ul- 
trarelativistic electron-positron plasma. It was shown that, in a dense ultrarel- 
ativistic electron-positron plasma electrostatic modes with wavelength greater 
than the plasma Debye length produce Langmiur Solitons that are spiky in na- 
ture. One essential feature of these solitons is that they can not form an energy 
flow toward smaller size, since they can not merge with each other. To merge, 
the soliton would have to give energy to sound waves, which do not exist in 
electron-positron plasma. In the one dimensional case the only process accruing 
in a soliton gas is nonlinear interaction with electron and positron, i. c. nonlin- 
ear Landau damping, which stops the solitons without changing its amplitude. 
The stopping length is of order of the soliton width, which is very small for a 
dense plasma. Thus the soliton will be stopped in a short time, of the order of 
their creation, these spiky short duration Langmiur solitons might berelated to 
the pulsar radiation and with its microstructures. 

In pi, V.I. Berezhiani, M.Y. El-Ashry and U.A. Mofiz are investigated the 
propagation of intense electromagnetic radiation in an admixture of unmagne- 
tized electron-positron-ion plasma analytically. It was shown that electromag- 
netic radiation of arbitrary amplitude in presence of heavy ions, in contrast to 
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the case of pure electron-positron plasma, may be localized with the generation 
of humped ambipolar potential in the plasma, i.e. the driving field intensity 
creates intense soliton in the plasma with the generation of double hump am- 
bipolar potentials. With the increase of the value of e {=n op /n oe ), they found 
that the tendency of converging a single hump soliton to a double hump one. 
In their investigation they neglected the ion dynamics. Consideration of ion 
dynamics, of course, may effect the localization phenomenon. This effect will 
be studied here. 

Now, the vector potential A for circularly polarized EM wave can be ex- 
pressed as: 

A = Aj_ + A z 

and 

A ± = a(z,t)e l ( fcz -^)+c.c. 
A z - a(z,t)e l ( fcz -^)+c. C . 

a(z,t) = a(z-Vt)(ic + iy) 
a(z, t) — d(z — Vt)z 

and a and d are real functions. 

To describe the admixture of plasma made up of electrons, positrons and 
ions, we use Maxwell equations, in which the fields are expressed in terms of the 
potentials, as described above. Accordingly, using Poissons equation, we may 
obtain the following equations for the potentials: 

Ax - c 2 V 2 Ax = 4tt C J - c—V4> (1) 

and: 

V 2 = -47T/9 (2) 

here, p and J are the charge and current densities, and given by: 

P = V" e a n a 



where: 



where a indicates the particle species a (=e, p, i for electrons, positrons, and 
ions, respectively), e a and n a are the charge and density of the corresponding 
particle a. We shall consider the case in which the admixture equilibrium state 
is characterized by n oe = n op + n i, where n oa , is the equilibrium density of the 
particle a. 

The relativistic equations of motion of different particles of unmagnetized 
plasma admixture is written as: 



d 



dt 



l<9Ax 



— P Q + m a c V7 Q = e a - V0 — (3) 



c dt 
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where 
and 



„2\-V2 / p2 \l/2 



where, m a is the rest mass of the particle a. 
The continuity equation for the particle a is: 

^ + V-(n a V a )=0 (4) 

We are looking for localized one-dimensional solution of this system of equa- 
tion for a circularly polarized EM wave, where the vector potential A can be 
expressed as: 

A = i(x + \y)A{C)e i{kz ~ ^ + c.c. 

here, £ = z — Vt. 

Analyzing the equation of motion and the wave equations into longitudinal 
(parallel to the direction of the wave propagation) and a transverse (perpendicu- 
lar plan) parts, then the transverse part of the equation of motion is immediately 
integrated to give: 

P a± = -^A ± (5) 

c 

where the constant of integration is set equal to zero, since the particles were 
assumed to be immobile at infinity where the field is zero. 

Meanwhile, the longitudinal part of the equation of motion takes the follow- 
ing form: 



dt +maC dz [1 



1/2 



OZ 



And, accordingly, equations( Q), and ( ||) may be rewritten as: 
d 2 1 d 2 \ . Aire 



(6) 



<)-.- ~ •-' m- ) A c ; " / ' V " 

+ n i V i -?i e V e ) (7) 

d 2 (j) 

-—- = -Aite(5n p + Sn, - 8n e ) (8) 
oz z 

Now, we consider only the electron; but everything stated below applies 
equally well, of course, to positron and ion. 
We have: 

d d d( d 
dt = d(dt = ~ V9 d( 
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and also 

d_ _ d_dC _ d_ 

dz d( dz d( 

then, equation( ^|) reduced to 



Now, we can integrate both sides with respect to (, yield: 



■2, /-, , P ez + P e± 



-v g P ez + mc\ 1 + 62 „ 2 e ^ = e0 + C 



where, C is the constant of integration, which could be determine from the 
boundary conditions. Here, P±,P ez — > Oand</> — > when £ — * 0, we get; C — 
mc 2 , and accordingly: 



(?) 



TOC 



Let: 



e0 



^ az — 



(9) 



c 



but for simplicity we will drop here the subscript e; i.e. <i> e = (E'and^'e = \& 
Equation ( p|) becomes 



-t/ 9 r e2 + v/i + r2 2 + * 2 = $ + i 

We note here that the began value of T ez is considered to be zero, then the 
above eaquation yields: 

r ez = iy<i> + i) 7g 2 - 

^y^+i) 2 - 1 ^ (io) 
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If we return back to the continuity equation; we get: 

dn e dv e 

U 71 

by integrating: 







-v g n e + n e v ez 



C 



here, v ez — ► when n e — > n oe , then, C = —v g n oe , and accordingly: 



U q n oe VT+T~fW 



but; 5n e = n e — n oe , hence 



5n P 



Uas/i + v ez + ^ - r e 



uwi + r ez +W -r e 



From equations ( and ( |ll|), we get: 



2 _ 1 + 

7 2 



here, we should note that: 



$ 

M 



where M is the mass of an ion. and that: 



m 
M 



(11) 



(12) 



By follows similar steps; for positron and ion; we can get the following: 

U g (l-$) 



5n p = n op j : 



(1 - $ 



|2 1+1' 2 



1 



(13) 



and 



§m = n oi % 



- 1 



(l-f*) 2 



(14) 
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The equation( pf) can be rewritten as: 



d 2 Aire 2 

■7— $ = T (5n e - 5n p - 5rii) 



Substituting equations( |T^), ( |T^) and ( |l4|) in equation( [l5| ): 



(15) 



(l + $) 2 



1+ip 2 

> a 



-ei 



Ug(l - $) 



(l-$) 2 --^ 

'a 



-C2 



(1 - - ~ 



M 

1 - ei - e 2 



But; 1 — ei — 62 = 1 
1 — ei — €2 — 0, hence: 



n oe —n v —n 



and n oe = n op + n oi ; we get, 



5?* = 



(l + $)2_ 1±*1 

fa 



-ei 



-C2 



t/,(l ^ *) 



( l_$)2_i±*: 



(l-f*) 2 - 



1+^* 2 



If we let ^ = ^ 



c 2 drj- 



■, we can write: 



<9 2 $ 
drj 2 



2r T 2 | 



-re/ 



(i + $) 



/ ( i + $ )2 _i±p 

ei(l-*) 



2 l+l' 2 



(i-$) 



1+^-* 2 



(16) 
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Equation ( ^) can be written as: 



dz 2 c? dt 2 ' J 1 



47re / n e n r , ran; i , 
+ — + tt~^ A A 



rat? V7e ' 7p ' M H 
By follows similar steps; the last equation can be written as: 



drf 




(17) 



Equation( |T^) and equation( [i~7]) are coupled together and forming a system 
of nonlinear second order deferential equations. It is rather difficult to obtain 
an analytical solution of them. However, in the following, we present numerical 
solutions of equation( [l(]) and equation( [l?]). Taking some typical parameters: 



0.12 
1 



M 1834 

Ug = VI -0.12 

and by using Taylor method, we have obtained the wakefield profiles at different 
values of ei and ti- The results are displayed in the following figures: 
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-30 



Fig2-a: The wakefield profiles for <I> at e\ = 0.4, e 2 = 0.6 
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